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Abstract. The issue of the X-ray to optical luminosity 
relationship [Lg — L^) is addressed for both optically and 
X-ray selected quasar samples. We have applied a general- 
ized regression algorithm for the case of samples involving 
censored data, with errors on both the dependent and the 
independent variable. Contrary to some previous results, 
we find that such relationship is consistent with being a 
linear one [L^ oc Lg). We argue that previous reports of 
non-linear relationships (i.e. oc with e < 1) are 
due to the neglect of the influence of the photometric er- 
rors, the precise knowledge of which strongly influences 
the reliability of the results. Further progresses in the de- 
termination of the Lo — Lx relationship can be achieved 
with ROSAT observations of the new generation of large 
bright quasar surveys. 
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^ 1. Introduction 

^ The study of the relationship between optical {Lg) and 
^ soft X-ray {L^) luminosities of quasars is needed to relate 
the quasar statistics (luminosity functions and evolution) 
in the two wavebands, and to better deflne the quasar's 
broad-band emission. So far the bulk of the available X- 
ray flux data has come from the Einstein mission. All 
comparative analyses of the evolution in the two bands 
have shown that the evolution rate is slower in the X-ray 
than in the optical (e.g. Maccacaro, Gioia &: Stocke 1984). 
Independently but consistently, a non-linear relationship 
between X-ray and optical luminosity of the type 



(1) 
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with e ~ 0.8, has been inferred by Kriss &: Canizares 1985, 
Avni &; Tananbaum 1986, Delia Ceca et al. 1992, Boyle 
et al. 1993. Recently Wilkes et al. (1994) have studied 
the Lo — Lx relationship in a sample of optically selected 
quasars observed by the IPC on board of Einstein, and 
have conflrmed the previous result by Avny &: Tananbaum 
(1986) of a non linear relationship between optical and X- 
ray luminosities. 

These questions have been reconsidered by Frances- 
chini et al. (1994), who showed that the evolution rate 
of X-ray selected AGNs could have been underestimated 
and actually be comparable to that of the quasars in the 
optical, with a linear scaling between the emission in the 
two bands. 

In this paper we explore in more detail the origin of 
the discrepancies on the Lo — L^ relationship, and discuss 
in particular the statistical uncertainties in the determi- 
nation of the regression parameters for both optical and 
X-ray selected samples. We flnally discuss the relevance of 
combining the ROSAT All Sky survey data with the new 
bright optical quasar samples, having a better photomet- 
ric quality than previous samples. We adopt throughout 
the values Ho = 50 Km s~^Mpc~^, Qo = for the Hubble 
and deceleration parameters. 

2. Generalized regression analyses 

Studying the relationship between the X-ray and optical 
luminosity of quasars, with the data presently available in 
the literature, is mainly a problem of flnding the correct 
regression algorithm and carefully determining the photo- 
metric errors. 

Different results on the relationship between two quan- 
tities are usually obtained from regression analyses as a 
consequence of which one is chosen as the dependent and 
which one as the independent variable. The simplest least- 
squares regression assumes that the independent variable 
is something given a priori and not subject to any obser- 



realistically take into account the fact that usually (as in 
our case) data points are affected by errors in both x and 
y coordinates. The slope of the best-fit straight line ob- 
tained in this way is systematically lower than the correct 
one. 

Two significantly different situations have to be dis- 
tinguished: X-ray selected samples or optically selected 
samples. 

In the first case, thanks to the potentially unlimited 
amount of observation time with ground-based telescopes, 
the optical follow-up is in practice complete, i.e. the op- 
tical data for the X-ray sources are not affected by cen- 
soring. Various statistical methods have been suggested 
to solve the problem in such a case (see e.g. Isobe et al., 
1990). In the following we will adopt a procedure of gen- 
eralized orthogonal regression, taking into account both 
measurements errors and intrinsic variances (i.e. scatter 
around the regression line unaccounted for by the obser- 
vational uncertainties) in the data (see also Fasano &: Vio 
1988). The problem is to estimate the parameters of the 
linear regression Y = a -\- bX, when also the dependent 
variable is subject to measurements errors. 

In terms of the first two empirical moments of the joint 
distribution F{x,y) of the data {x} and {y\ we have: 
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Here x, y are the sample means, cov[x, y) is the covariance 
and (7^ is the variance of the distribution of the x variable: 

cov{x, y) =< {x - x){y - y) >, 

al =< {x -xy>. 

If also the dependent variable is subject to experimen- 
tal errors, the estimator of the slope in eq. (2) is biased 
toward a lower value. Let's define and Cy the errors in 
the X and y axes, and (7e^ and (Te^ their standard devia- 
tions, with (7e^ ^ (Te^ but Constant for all the data. In this 
case (7^ = o"^ + (7g where al is the "true" variance of the 
distribution of x and a^^ is the contribution due to the 
measurement errors. Since the expected value oi cov[x,y) 
is not affected by the measurement errors of x, the larger 
is (7g^ the lower will be the expected value of b. If b is 
biased, according to eq. (3) also the estimator a will be 
biased. 

Assuming gaussian errors, we can associate to each ob- 
served point [xi,yi) a hivariate reormaZ probability density 
function pi[x,y), whose characteristic variances along the 
axes are a^^ and a^^. In this case, the regions around each 
observed point, where the "true" point can be found with 
some given probabilities (confidence regions), are bounded 
by the families of ellipses 



(4) 



at varying fc,- . The point of the true straight line that has 
most likely generated a given observed point is that tan- 
gential to the family (4). The tangency condition is found 
to be 



= iVi -a-bxifl{al 



a. 



According to the Maximum Likelihood principle, the most 
likely straight line is obtained by maximizing the product 
of the probabilities pi[x,y) computed at the tangential 
points, that is by minimizing the sum 



(5) 



The minimization of quantity (5) with respect to a and b 
provides the estimator in eq. (3) and the equation 
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It is possible to show that the same results are obtain- 
able by following a Least Squares approach (see Press &: 
Teukolsky, 1992). 

In the case of optically selected samples, because of the 
limits in the duration of X-ray missions, in addition to the 
uncertainties due to the errors on both axes, the presence 
of upper limits to the X-ray flux has to be taken into ac- 
count. Avni et al. (1980) and Avni e Tananbaum (1986) 
have developed a regression method which deals with up- 
per limits, that is based on the principles of the "Detec- 
tions and Bounds" (DB) method, but does not consider 
the influence of errors on both axes. For the most general- 
ized cases, in which both upper limits and errors in the x 
and y axes are present, we have used a simple extension of 
the Schmitt's technique (Schmitt, 1985). Note that other 
linear regression methods are available (for a review see 
Isobe, Feigelson &: Nelson 1986), but their extension to 
the case of errors in both axes is not straightforward. 

In its original formulation, the Schmitt's method pro- 
vides an estimate of the parameters of the linear regres- 
sion y = a + 6Jl , computing the estimators (2) and (3) 
by taking into account the presence of upper limits in the 
computation of x, y, coi)(a;,i/), a^.. It is then straightfor- 
ward to take into account the presence of errors on both 
axes by substituting the estimator (2) with estimator (6), 
and by computing the mean values x, y, the covariance 
and the variance of distribution of the x variable, exactly 
as explained by Schmitt (1985). 



AGNs 

The Lq — Lx relationship for X-ray selected AGNs has been 
extensively investigated by Franceschini et al. (1994). We 
outline here the main results. The sample used is the com- 
bination of the EMSS (Stocke et al. 1991) with the deep 
ROSAT observations from Boyle at al. 1993, resulting in 
the largest sample available at present. 

A generalized regression analysis allowing for errors 
and intrinsic variances in both Mb and has been 
performed as explained in Section 2. We have assumed 
that the Einstein X-ray fluxes have uncertainties of 30%, 
a value somewhat larger than reported by Gioia et al. 
(1990), but realistically closer to account for various addi- 
tional uncertainties, in particular those due to the effect of 
a variable galactic low energy absorption. An error of 0.5 
magnitudes on Mb has also been assumed. As explained 
below, this takes into account: 1) that the EMSS optical 
photometry is a collection of CCD data and Schmidt plate 
POSS measurements (with an average error of 0.4 magni- 
tudes), 2) that the quasar variability adds an uncertainty 
of ~ 0.2 magnitudes, 3) and that an error of 0.15 magni- 
tudes due to intrinsic spread on the optical spectral index 
should be taken into account. 

The effect of the variability is difficult to address, as 
it significantly depends on the quasar absolute optical lu- 
minosity, and increases with decreasing redshift. From the 
analysis of the structure function (Hook et al. 1994) the 
typical amplitude of the quasar variability increases with 
the time-base in the quasar rest frame according to the 
formula: 

< |Am| >= [0.196 + 0.033(Mb + 25.7)](1 - e'^l'^-^^fl'^ , 

(7) 

where r is expressed in years. Even if we take into ac- 
count the redshift correction, it results that on average 
the optical photometry of the Einstein AGNs has been 
carried out some r > 0.38(1 + z) years later or before 
the X-ray measurements (at a typical redshift z ~ 0.5 for 
the EMSS), and we should expect an average error due to 
variability of 0.2 magnitudes for quasar having a typical 
absolute magnitude Mb ~ —26 (0.1 for Mb ~ —30 and 
0.3 for Mb 23). 

The dispersion of the optical spectral index a [Fu oc 
is typically 0.25-0.30 (e.g. Sargent, Steidel &: Bok- 
senberg, 1989), this figure originating an error of ~ 0.15 
magnitudes at an average redshift ~ 1, in agreement with 
what observed by La Franca, Cristiani &: Barbieri (1992). 

We have excluded from our analysis objects having 
Lx{o.3-3.z keV) < lO^^ergfs"^ and z < 0.2 in order to 
avoid the presence of spurious trends introduced by the 
contribution of the low-redshift, low-luminosity objects in 
which the emission from the host galaxy becomes non- 
negligible. It results that the functional dependence be- 
tween the X-ray and optical luminosity for AGN at z > 0.2 
and Lx > 10^^ is basically a linear one, with no indication 



tical to X-ray luminosity ratio. The best-fit regression for 
268 AGNs is 

Mb = -2.5(±0.1) logL(^o.3-3.5 keV) + 87.9(±4.9), (8) 



(9) 



(quoted uncertainties are bootstrap errors at 95% confi- 
dence). This corresponds to: 



log{Lo/Lx) _ log{Lx/Lo) 
log{vxlvo) 2.605 



1.32, 



(10) 



L(2 keV) and logL^^^^^ being expressed in units of 

erg/ s/ Hz. An essentially linear relationship for medium- 
to-high luminosity objects is confirmed by a detailed anal- 
ysis of the residual Mb distribution. In particular, by 
adopting our best-regression solution, a Gaussian model 
with (7(AMb) = 0.85 fitted to the binned distribution of 
the residuals in eq. (8) provides a good fit (with a total 
= 21 for 22 independent bins). 



4. The Lo 
AGNs 



Lx relationship for optically selected 



The relationship between optical and X-ray emission for 
optically selected AGNs has been studied by various 
groups using data from the Einstein mission. Avni and 
Tananbaum (1986) using a sample with 94 X-ray detec- 
tions and 60 upper limits found that the relationship be- 
tween X-ray and optical luminosity was non-linear, with 
a best-fit value for the exponent in eq. 1 of e ~ 0.8. More 
recently and for a richer sample (179 detections and 164 
upper limits), Wilkes et al. (1994) found a best-fit value 
for the exponent as small as e ~ 0.7. 

We have verified that these dependences cannot apply 
to X-ray selected samples. Regressions lines with e < 1 en- 
tail residual distributions both markedly correlating with 
the independent variable and characterized by a signif- 
icant skewness. If we adopt in particular e = 0.7 in eq. (1), 
we would obtain a distribution of the residuals which is 
clearly a function of the X-ray luminosity and to which a 
Gaussian model (even with a standard deviation increased 
to (7(M) = 0.95) would provide a very poor fit. 

In order to better understand the inconsistencies be- 
tween results of regression analyses of X-ray selected sam- 
ples and those obtained for optically selected samples, we 
have analyzed the data of Wilkes et al. (1994) applying 
the method for DB data with errors on both axes as ex- 
plained in Section 2. As discussed there, the results de- 
pend on the assumed errors of the measurements on the 
two axes. Unfortunately, the uncertainties on and Mb 
are not known, but it seems reasonable to consider them 
comparable in the logLg — logL^ plane: an error of 30% 
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Fig. 1. Correlation between optical luminosity at 2500 A and 
X-ray luminosity (2 keV) for optically selected samples. The 
filled circles indicate detections, the open circles indicate upper 
limits, (a) The sample of Wilkes et al. (1994): the continuous 
line is our best fit using the DB method, the dotted lines are 
the two extreme slopes in the case in which errors are assumed 
in turn only on the x or on the y axis; the dashed line is the 
best fit obtained by Wilkes et al. (1994). (b) Simulation of 
ROSAT observations (see text). 46 detections (filled circles) 
and 80 upper limits (open circles). The range in luminosities is 
smaller than in case a) because of the exclusion of the objects 
having Mb > -23. 



in the fluxes as used for the Einstein data in the previ- 
ous section corresponds to 0.75 magnitudes. Under this 
assumption we have obtained the relation 

logL^2 keV) = 1.02(±0.04) logL^^^^^^^ - 4.5(±1.1), (11) 

(quoted uncertainties are bootstrap errors at 95% confl- 
dence) in agreement with the linear dependence between 
the two quantities found in Section 3 for the X-ray se- 



errors are assumed in the optical fluxes (this correspond 
to the result found by Wilkes et al. 1994) and 1.35 if no 
errors are assumed in the X-ray fluxes (see Figure la). In 
this case, following Wilkes et al. (1994), also the objects 
having ix(o.3-3.5 keV) < 10^^ erg have been included 
in the regression analysis. 

As discussed by Cheng et al. (1984), X-ray flux limited 
observations will preferentially select objects with larger 
values of L^/Lo, i.e. with smaller values of agx- This ex- 
plains the different constant terms in the Lg — rela- 
tionships for X-ray selected (eq. [9]) and optically-selected 
quasars (eq. [11]). Let's indicate with a.ox,x and cxox, Optical 
the spectral indices computed using X-ray and optically 
selected samples respectively, and with (Ta^x the standard 
deviation of the "observed" distribution oiaox,x- Accord- 
ing to Cheng et al. (1984) and assuming a Gaussian dis- 
tribution of the residuals {cmb = 0.85, or (Ta^x = 0.153), 
we obtain 



aox,Opticai = oiox,x + 90-2 ^ = 1.32 + 9(0.153)^ = 1.53, 

(12) 

in good agreement with cxox, Optical = 1-49 obtained from 
eq. (11) at logLg = 31. 

5. Discussion and Conclusions 

Our analysis shows that a linear relationship between opti- 
cal and X-ray luminosity seems to be required by the avail- 
able data from X-ray selected samples, and consistent with 
those on optical samples. Previous results favoring non- 
linear relationships may be due to the underestimation of 
the errors in the optical fluxes. Additional uncertainties 
are introduced by the large number of upper limits to the 
X-ray fluxes in the optical samples, by emission and/or 
extinction of the optical light from the host-galaxy in low- 
luminosity objects, and by the sparse nature of optical 
samples. From the above discussion it results clear that in 
order to unequivocally address the Lg — Lx relationship it 
is necessary to analyze statistically-well-deflned large data 
samples not only with an high level of photometric accu- 
racy, but also with an accurate knowledge of the errors 
affecting the measurements on both the X-ray and opti- 
cal fluxes. If the time distance between the measurements 
in the two bands is larger than few years, the variability 
adds errors of ~ 0.1 — 0.3 magnitudes. The achievement of 
accurate optical photometry provides little improvement 
if these effects are ignored. 

Good opportunities in this sense are provided by recent 
advances: on the one hand the advent of the ROSAT satel- 
lite is providing new X-ray data, on the other hand three 
new bright optical quasar surveys are becoming available: 
the Edinburgh quasar survey (ROE, Goldschmidt et al. 
1992); the Large Bright Quasar Survey ( LBQS, Hewett, 
Foltz &: Chaffee 1993 and references therein); the Homoge- 



1995). 

We have analyzed the level of accuracy in the determi- 
nation of the Lo — Lx relationship that could be achieved 
from the ROSAT "All Sky Survey" data of quasars from 
two of these surveys: the HBQS and ROE samples, having 
almost the same color selection criterion, and the highest 
photometric quality. As for the analysis of the X-ray se- 
lected quasars in Section 3, we have excluded the fainter 
objects having Mb > —23, in order to avoid spurious 
trends introduced by the contribution of the light from 
the hosting galaxy, and the incompleteness which affects 
color selections for the intrinsically fainter quasars. 

The HBQS and ROE surveys give a complete sample of 
about 120 quasars over an area of more than 800 deg^ with 
15 < 5 < 17.5 (the inclusion of the LBQS would bring 
the sample size to roughly 200 quasars in this magnitude 
range). Particular care has been dedicated to the photo- 
metric accuracy (e.g. the correction of dis-uniformities in 
the response of the photographic material within a given 
Schmidt plate) reaching absolute errors less than 0.1 mag. 
To compute the total error on the optical fluxes which will 
be compared with the X-ray fluxes, we should also take 
into account the effect of variability (eq. 7), and of the 
dispersion of the optical spectral index as done in Section 
3. For the HBQS there is a typical difference of about two 
years between the observations in the optical and ROSAT 
bands. The ROE survey has larger time differences. As- 
suming an average redshift of z ~ 1, according to eq. (7) 
we can estimate an error due to variability ~ 0.2 mag- 
nitudes for both surveys. As discussed in Section 3, the 
dispersion of the optical spectral indices causes an error 
of ~ 0.15 magnitudes at an average redshift z ~ 1. Then 
we expect a typical total uncertainty on the optical mag- 
nitudes (when reported at the time of the X-ray observa- 
tions) of (7b ~ 0.25 mag. The advantage in using samples 
such as the HBQS and ROE is to allow a well deflned sta- 
tistical deflnition and treatment of the photometric errors 
which represents a considerable improvement with respect 
to previous analyses which were based on incomplete sam- 
ples from published inhomogeneous catalogs of quasars. 

To estimate the errors on the X-ray measurements, 
a linear regression between the S/N ratio and the X-ray 
-fx(o.5-2 keV) ROSAT fluxes has been computed from a 
sparse sample of more than 200 quasars selected from the 
Veron &: Veron (1993) catalogue and observed by the "All 
Sky Survey". It results: 

log{ — ) = 0.53(±0.02)/o3F^ + 7.09(±0.08) (13) 

(Molendi &: Doublier, private communication). We have 
consequently used an "All Sky Survey" 3(7 average de- 
tection limit of Fj;(o. 5-2.5 fcei^) = 3 10~^^ erg cm~^s~^, 
and have simulated the almost linear relationship be- 
tween Lx and Lo as in eq. (11) for the quasars of the 
HBQS and ROE samples. Eq. (13) gives an average error 



HBQS and ROE surveys. 

A plot of Lx versus Lo for simulated ROSAT observa- 
tions of the HBQS and ROE quasars appears in Fig. lb. It 
results that roughly one third of the quasars are expected 
to have a detection at a 3(7 confldence level. Assuming 
(7b = 0.25 and (yiogF^ = 0.17 constant for all the data, our 
method of DB allowing for errors on both axes reproduces 
the assumed value for the slope (eq. 11): h = 1.02 ± 0.03. 
In the extreme cases in which the errors are assumed only 
either on the X-ray or on the optical fluxes, the slope 
value varies from 0.95 to 1.09. This is a quite narrower 
range than the interval 0.76 < e < 1.35 previously ob- 
tained from the data collection by Wilkes et al. (1994), 
and shows the potential accuracy that could be achieved 
from this kind of observations. 

Added to these uncertainties we should expect an in- 
trinsic spread in the Lo — L^ relationship, which at present 
is completely hidden by the the photometric errors in the 
data collected by Wilkes et al. (1994). Exploitation of the 
ROSAT "All Sky Survey" data for quasars in the HBQS 
and ROE samples will allow to test this idea. 

The symmetric approach to the problem, based on op- 
tical identiflcations of ROSAT X-ray selected samples, is 
underway (Danziger et al., 1990, and Barcons et al., pri- 
vate communication). As demonstrated above, also in this 
case the reliability of the results will be strongly depen- 
dent on the precise knowledge of the photometric errors. 
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